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The i 2 -moment of the twist-3 polarized parton distribution gi(x) is related to the transverse force 
acting on the active quark in deep-inelastic scattering off a transversely polarized nucleon imme- 
diately after it has absorbed the virtual photon. Lattice calculations of the corresponding matrix 
element as well as experimental measurements of g2{x) are used to estimate sign and magnitude of 
this force. Similarly, the i 2 -moment of the chirally odd twist-3 unpolarized parton distribution e(x) 
can be related to the transverse force experienced by a transversely polarized quark ejected from a 
transversely polarized nucleon. 

PACS numbers: 



I. INTRODUCTION 

In the analysis of high-energy scattering processes, leading twist effects are often both easier to isolate (e.g. by 
increasing Q 2 until higher twist effects are negligible) and have a more direct physical interpretation than contributions 
from higher twist operators which are often intermingled with -q^ corrections to leading twist operators. A notable 
exception is is the twist-3 polarized structure function g2(x,Q 2 ), which can be cleanly separated from its twist-2 
counterpart gi(x, Q 2 ) when scattering longitudinally polarized leptons from a transversely polarized target. While in 
general, the contribution from g^ to the polarized cross section is suppressed by powers of , at a target polarization 
angle of 90 ° the leading contribution vanishes and the contribution from g\ is exposed without kinematical suppression. 
This property of the polarized DIS cross section thus allows a clean extraction of higher twist matrix elements, without 
the need for fitting and subtracting a leading twist effect, and thus makes polarized deep-inelastic scattering (DIS) a 
rare opportunity for studying higher twist effects (for an overview, see Ref. [l[). 

As g2(x,Q 2 ) involves higher twist, it does not have a parton interpretation as a single particle density. Indeed, 
the twist-3 part of gi is related to quark-gluon correlations whose intuitive interpretation may not be immediately 
clear. Since g2(x,Q 2 ) is related to (electromagnetic) polarizabilities at low Q 2 , these twist-3 matrix elements have 
been called color polarizabilities in the literature 0. However, at high Q 2 , the twist-3 piece of g2(x, Q 2 ) is described 
by a local correlator and the physical interpretation as a polarizability no longer applies. Indeed, while nucleons need 
to be polarized in order to study g%(x, Q 2 ), the nucleons are not distorted, but only 'spin-alligned'. The quark-gluon 
correlations embodied in the twist-3 part of g2{x,Q 2 ) are then obtained as a matrix elements of a certain operator 
in a spin-alligned, but undcformed, nucleon. This is very different from the usual use of the term 'polarizability' as 
the tendency of a charge or magnetization distribution to be distorted from its normal shape by an external field. Of 
course, one could broaden the notion of 'polarizability' to encompass matrix elements that are only non-zero when the 
nucleon is polarized, but within such a broadened definition, other spin-dependent observables, such as the polarized 
parton distribution Aq(x) or even the magnetic moment of the nucleon, would then also become 'polarizabilities' in 
the broader sense. 

The main purpose of this paper is to explore an alternative physical interpretation of these particular twist-3 
matrix elements as a force. First we summarize the connection between the x 2 moment of 32(^7 Q 2 ) and quark-gluon 
correlations. After discussing the connection between these correlations and the transverse force on the active quark in 
DIS, we then estimate sign and magnitude of that force based on DIS data, lattice calculations and heuristic pictures. 

II. x 2 MOMENTS AND QUARK-GLUON CORRELATIONS 

The chirally even spin-dependent twist-3 parton distribution g2{x) — gr(x) — gi(x) is defined as 

— e lXx (PS\mn 5 ^n)\ Q 2 \PS) = 2 [ 9l (x, Q 2 )p^(S ■ n) + g T (x, Q 2 )S^ + M 2 g 3 (x, Q 2 )n»(S ■ n)] . (1) 

where p^ and are light-like vectors along the — and + light-cone direction with p ■ n — 1. Using the equations of 
motion 52 (x) can be expressed as a sum of a piece that is entirely determined in terms of g\ (x) plus an interaction 
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dependent twist-3 part that involves quark gluon correlations Q 

92(x) = g 2 vw (x)+g 2 (x) 



(2) 



g? w (x) = -gi(x)+ / -Z- gi (y). 

j x y 

Here we have neglected m q for simplicity. For example, the x 2 moment yields [3, [Bj] 



/ 



dxx 2 g 2 (x) = 



with 



g (P, S \q(0)G +y (0)j+q(Q) \P,S)= 2MP+P+S x d 2 



(3) 



(4) 



Note that conventions exist where ^ appears on the r.h.s. of (J3]) in which case there is no factor 2 on the r.h.s. of 
©■ 

In the limit where Q 2 is so low that the virtual photon wavelength is larger than the nucleon size, the electro- 
magnetic field associated with the two virtual photons appearing in the forward Compton amplitude corresponding to 
the structure function is nearly homogenous accross the nucleon and the spin-dependent structure function g 2 {x, Q 2 ) 
can be related to spin-dependent polarizabilities. In contradistinction, in the Bjorkcn limit, the matrix elements 
describing the moments of g 2 (x, Q 2 ) are given by local correlation functions, such as (g]). Nevertheless, because of the 
abovementioned low Q 2 interpretation of g 2 , the local matrix elements appearing in ((4]) 



XE 2M 2 S = (P, S\ q^a X gEq \P, S) 



X b2M 2 S= (P,S\qigBq\P,S) 



(5) 



where 



(Xe ~ 2 xm) , 



(6) 



(note that ^/2G +V = B x - E») are sometimes called color electric and magnetic polarizabilities [2]. In the following 
we will discuss why, at high Q 2 , a better interpretation for these matrix elements is that of a 'force'. 

In electro-magnetism, the y-component of the Lorentz force F y acting on a particle with charge e moving, with 
(nearly) the speed of light along the — z direction, v as (0, 0,-1), reads 



py — 



E + vxb] = e(E v - B x ) = -eV2F +y , 



(7) 



which involves the same linear combination of Lorentz components that also appears in the gluon field strength tensor 
in j?]) . This simple observation already suggests a connection between d 2 and the color Lorentz force on a quark that 
moves (in a DIS experiment) with v ~ (0, 0, — 1). 

In order to explore this connection further we compare the matrix element defining d 2 with that describing the 
average transverse momentum of quarks in semi- inclusive DIS (SIDIS) [□]. The average intrinsic transverse momentum 
of quarks bound in a nucleon vanishes and therefore any net transverse momentum of quarks in a SIDIS experiment 
must come from the final state interactions (FSI) [?J • The average transverse momentum of the ejected quark (also 
averaged over the momentum fraction x carried by the active quark in order to render the matrix element local in 
the position of the quark field operator) in a SIDIS experiment can thus be represented by the matrix element @ 



(k y ± ) 



1 



2i^ 



P,S 



q(0) dx-gG +y (x+ = 0,x-)j+q(0) 



P,S 



(8) 



where Wilson-line gauge links along x~ are implicitly understood, but not written out explicitly. One way to derive 
this expression is to start from gauge invariantly defined quark momentum distributions with Wilson line gauge links 
extending to ligh-cone infinity. The integral over the gauge field is then obtained by acting with the transverse 
derivative (when averaging over fc^ ) on the gauge field appearing in the Wilson line 0, [Ic| . 

The matrix element appearing in ((SJ) thus has a simple physical interpretation as the transverse impulse obtained by 
intergrating the color Lorentz force along the trajectory of the active quark — which is an almost light-like trajectory 
along the —z direction, with z = —t. 
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In order to make the correspondence more explicit, we now rewrite Eq. ((8|) as an integral over time 

(kl) = -2p+ ( p > s \m J dtG+y(t } z = -t) 7 +q(0) |P, 5) (9) 

in which the physical interpretation of —jpr (P, S\ q(0)G +v (t, z = —t)~/ + q(0) |P, S) as being the (ensemble averaged) 
transverse force acting on the struck quark at time t after being struck by the virtual photon becomes more apparent. 
In particular, 

^(0) = (P,S\q(0)G +y (0) 7 + q(0)\P,S) (10) 

= -V2MP + S x d 2 = -M 2 d 2 , 

where the last equality holds only in the rest frame (p + = -^M) and for S x — 1, can be interpreted as the averaged 
transverse force acting on the active quark in the instant right after it has been struck by the virtual photon. 

Although the identification of {p\<Ti + G +v q\p) as a color Lorentz force may be intuitively evident after the above 
discussion, it is also instructive to provide a more formal justification. For this purpose, we consider the time 
dependence of the transverse momentum of the 'good' component of the quark fields (the component relevant for DIS 
in the Bjorken limit) q+ = \^~~l + q 

2p + j t {p y ) = j t (PS\ g 7 + ( P y - gAy) q \PS) = -±=j f (PS\ q\ { P y gA«) q+ \PS) (11) 

= 2 P + (ps\ { v y - 9 Ay) q + qj+ ( p y - g Ay) q - g 7 + <?i^] \ps) . 

Using the QCD equations of motion 

q= (igA°+ 7 °y-D)q, (12) 

where —iD^ = p^ — gA 11 , yields 

2p + ^(p y > = {PS\q 1 + g(Gy a + Gy*)q\PS) + '{PS\q 1 + 1 - 1 *D l Diq\PS)' (13) 
= y/2 (PS\ qi + gG y+ q \PS) + ' (PS\ qy + 'Y~i i D i B i '<1 \PS)' , (14) 

where ' (PS\ qj + '7' D l D 3 q \PS)' stands symbolically for all terms that involve a product of 7 + 7~ as well as a 7" 1 
and that also involve only transverse derivatives D l . 

Now it is important to keep in mind that we are not interested in the average force on the 'original' quark fields 
(before the quark is struck by the virtual photon) , but after absorbing the virtual photon and moving with (nearly) 
the speed of light in the — z direction. In this limit, the first term on the r.h.s. of (114|) dominates, as it contains the 
largest number of '+' Lorentz indices. Dropping the other terms yields (|10[) . 

A measurement of the x 2 -moment f% of the twist-4 distribution g^(x) allows determination of the expectation value 
of a different linear combination of Lorentz/Dirac components of the quark-gluon correlator appearing in ((J) [Tlj 

f 2 M 2 S» = \ (p, S\ qg&»w \p, S) . (15) 

Using rotational invariance, to relate various Lorentz components one thus finds a linear combination of the matrix 
elements of electric and magnetic quark-gluon correlators ([5]) 

h = XE- Xm, (16) 

that differs from that in ([6]). In combination with ([4]) this allows a decomposition of the force into electric and 
magnetic components F v = P| + F\ I: using 



n(0) = --rXE f v b (0) = -—xb (17) 



M 2 11 / \ M 2 

—XE F y B (0) = -^ 

for a target nucleon polarized in the +x direction, where 0, [l2| 



Xe = I (2d 2 + f 2 ) xm = I (4d 2 - /a) • (18) 
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FIG. 1: Distribution of the j + density for u and d quarks in the _L plane (xbj = 0.3 is fixed) for a nucleon that is polarized in 
the x direction in the model from Ref. [151 ]. For other values of x the distortion looks similar. 



A relation similar to (p~0|) can be derived for the x 2 moment of the twist-3 scalar PDF e(x). For its interaction 
dependent twist-3 part e(x) one finds for an unpolarized target fl3l ] 

2 

4MP+P+e 2 = ^2g(p\qa +l G +l q\P) , (19) 

i=l 

where e 2 = Jq dxx 2 e(x). The matrix element on the r.h.s. of Eq. (fTi?)) can be related to the average transverse 
force acting on a transversely polarized quark in an unpolarized target right after being struck by the virtual photon. 
Indeed, for the average transverse momentum in the +y direction, for a quark polarized in the +x direction (kj_ 
moment of the Boer- Mulders function [l4[ integrated also over x) , one finds 

1 f°° 

(k y ) = ^tJ q (Pi <j(0)<T +y G +y (x-)q(0) \p) . (20) 

A comparison with Eq. (fl"9|) shows that the average transverse force at t = (right after being struck) on a quark 
polarized in the +x direction reads 

1 1 M 2 

F y (0) = ^-j^g (p| qa+yG+yq \p) = —MP+S*e 2 = — e 2 , (21) 

where the last identify holds only in the rest frame of the target nucleon and for S x = 1 . In the physical interpretation 
of (|21|) it is important to keep in mind that, for a given flavor, the number of quarks on which the force in (121|) is 
only half that in (|10|) as only half the quarks in an unpolarized nucleon will be polarized in the +x direction. 



III. HEURISTIC PICTURES AND NUMERICAL STUDIES 



When the target nucleon is transversely polarized, e.g. in the +x direction the axial symmetry in the transverse 
plane is broken. In particular, the quark distribution (more precisel y th e distribution of the 7 + -density that dominates 
in DIS in the Bjorken limit) in the transverse plane is deformed The average deformations can be related to 

the contribution from each quark flavor to the anomalous magnetic moment of the nucleon and was predicted to be 
quite substantial [l5| and has also been observed in lattice QCD [l6|]. Given the fact that, for a nucleon polarized 
in the +x direction the 7 + -distribution for u (d) is shifted towards the ±y direction suggests that these quarks also 
'feel' a nonzero color-electric force pointing on average in the =Fy direction, i.e. one would expect that d 2 is positive 
(negative) for u (d) quarks. This is also consistent with a negative (positive) sign for the Sivers on the proton function 
as observed by the HERMES collaboration [TtJ and the vanishing Sivers function for deuterium in the COMPASS 
experiment [l8(. In fact, in the large Nc limit, one would expect that d 2 for u (d) quarks are equal and opposite. 
Note that while d 2 for u (d) quarks being exactly equal and opposite would imply the same for protons (neutrons), 
any deviation from being exactly equal and opposite is enhanced for proton (neutron) since there is a significant 
cancellation between the two quark flavors in the nucleon. 

If all spectators in the FSI were to 'pull' in the same direction, the force on the active quark would be of the order 
of the QCD string tension a w (45QMeV) 2 , which would translate into a value d 2 ~ 0.2. However, it is more natural 
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to expect a significant cancellation between forces from spectators pulling the active quark in different directions, 
the actual value of d^ is probably about one order of magnitude smaller, i.e. d-i ~ 0.02 appears to be more natural. 
Instanton based models have suggested an even smaller value (l9j |. 

Heuristic arguments/lattice calculations [2(| also suggest that the deformation of (the 7 + -distribution for) 
transversely polarized quarks in an unpolarized nucleon is more significant than that of unpolarized quarks in a 
transversely polarized nucleon. When applied to the final state interactions, this observation suggests \ei\ > \ds\ (the 
fact that in an unpolarized nucleon only half the quarks are polarized in the -fi-direction is compensated by the 
factor | in (|2Tj) . 

Lattice calculations of the twist-3 matrix element yield [2l[ 

4 u) = 0.020 ± 0.024 4 d) = -0.011 ± 0.010 (22) 

renormalized at a scale of Q 2 = 5 GeV 2 for the smallest lattice spacing in Ref. [2l| . Note that we have multiplied the 
numerical results from [2l| by a factor of 2 to account for the different convention for di being used. 
Here the identity M 2 ss 5GeV/fm is useful to better visualize the magnitude of the force. 

F {u) = -25 ± 30MeV/fm F [d) = 14 ± 13MeV/fm. (23) 

In the chromodynamic lensing picture, one would have expected that Fr u \ and are of about the same magnitude 
and with opposite sign. The same holds in the large Nc limit. A vanishing Sivers effect for an isoscalar target would 
be more consistent with equal and opposite average forces. However, since the error bars for di include only statistical 
errors, the lattice result may not be inconsistent with ~ — d% ■ 

The average transverse momentum from the Sivers effect is obtained by integrating the transverse force to infinity 
(along a light-like trajectory) (k v ) = J °° dtF y (t) This motivates us to define an 'effective range' 

Reff = FWr ( } 

Note that R e f / depends on how rapidly the correlations fall off along a light-like direction and it may thus be larger 
than the (spacelike) radius of a hadron. Of cource, unless the functional form of the integrand is known, R e ff cannot 
really tell us about the range of the FSI, but if the integrand does not oscillate 

Fits of the Sivers function to SIDIS data yield [H| one finds about | (k v ) | ~ 100 MeV [II]. Together with the (average) 
value for |g?2| from the littice this translates into an effective range R e ff of several fen. It would be interesting to 
compare R e ff for different quark flavors and as a function of Q 2 , but this requires more precise values for d2 as well 
as the Sivers function. 

Note that a complementary approach to the effective range was chosen in Ref. [23], where the twist-3 matrix 
element appearing in Eq. (fTOjl was, due to the lack of lattice QCD results, estimated using QCD sum rule techniques. 
Moreover, the 'range' was taken as a model input parameter to estimate the magnitude of the Sivers function. 

The im pac t parameter distribution for quarks polarized in the +x direction was found to be shifted in the +y 
direction jl6l . [20 , [24| . Applying the chromodynamic lensing model implies a force in the negative —y direction for 
these quarks and one thus expects ei < for both u and d quarks. Magnitude: since kj_ > k, expect odd force larger 
than even force and thus | e2 1 > | 1 - 

It would be interesting to study not only whether the effective range is flavor dependent, but also whether there is 
a difference between the chirally even and odd cases. It would also be very interesting to learn more about the time 
dependence of the FSI by calculating matrix elements of qj + (Z? + G +± ) q, or even higher derivatives, in lattice QCD. 
Knowledge of not only the value of the integrand at the origin, but also its slope and curvature at that point, would 
be very useful for estimating the integral in Eq. ([H]). 



IV. DISCUSSION 



The quark-gluon correlations in the a; 2 -moment c?2 of the twist-3 polarized PDF can be identified with the 
transverse component of the color-Lorentz force acting on the struck quark in the instant after absorbing the virtual 
photon. The direction of the the force for u and d quarks can be understood in terms of the transverse deformation of 
parton distributions for a transversely polarized target. In combination with a measurement of the x 2 moment of the 
twist-4 polarized PDF 173 one can even decompose this force into color-electric and magnetic components. Although 
still quite uncertain, first eperimental/lattice results suggest values around 25 — 50MeV/fm for the net force. This 
should be compared with the net transverse momentum due to the Sivers effect which is on the order of lOOMeV. 
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The x 2 moment e-i of the chirally odd twist-3 (scalar) PDF e(x) can be related to the transverse force acting on 
transversely polarized quarks in an unpolarized target. Therefore, e 2 is to the Boer-Mulders function hj^, what is g?2 
to the Sivers function f^ T . 
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